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Vortex flow remains laminar up to large (~ 10^) Reynolds numbers Rca in a cylinder filled with 
^He-B. This is inferred from NMR measurements and numerical vortex filament calculations where 
we study the spin up and spin down responses of the superfluid component, after a sudden change 
in rotation velocity. In normal fluids and in superfluid *He these responses are turbulent. In ^He-B 
the vortex core radius is much larger which reduces both surface pinning and vortex reconnections, 
the phenomena, which enhance vortex bending and the creation of turbulent tangles. Thus the 
origin for the greater stability of vortex flow in ^He-B is a quantum phenomenon. Only large flow 
perturbations are found to make the responses turbulent, such as the walls of a cubic container or 
the presence of invasive measuring probes inside the container. 

PACS numbers: 67.30.hb, 47.15.ki, 02.70.Pt, 67.30.he 
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For centuries the transition to turbulence has been one 
of the enigmatic unsolved problems of classical hydrody- 
namics. It is commonly accepted that laminar flow be- 
comes unstable with increasing Reynolds number Rei, 
above some case-dependent threshold Reth ^ 1...10^. 
This instability results either in secondary vortical flow 
or (depending on the flow geometry) in intensive tur- 
bulent motions with many degrees of freedom. In spe- 
cial geometries, particularly in a circular pipe, laminar 
flow is assumed asymptotically stable with respect to in- 
finitesimal perturbations at any Rcj^ ^, but the critical 
amplitude of perturbation, which triggers turbulence, de- 
creases rapidly as ex l/Re^ [2]. In practice, normal fluid 
flows are thus turbulent for Re^, ^ Reth- 

What can be said about the transition to turbulence 
in the flow of quantized vortices in superfluids? When 
the temperature T ^ 0, dissipation from conventional 
mechanisms approaches zero and laminar vortex flow is 
believed to become unstable. This is the current impres- 
sion about 3-dimensional vortex flow in superfluid ^He. 
Here we show that in superfluid ^He-B laminar vortex 
flow can be stable down to temperatures of 0.2 Tc in an 
axially symmetric situation and, to make it turbulent, a 
strong perturbation is needed. 

Laminar vortex flow:-The Reynolds number Re^, = 
VL/v provides an estimate of the ratio of the nonlinear 
(^ V'^ /L) and dissipative (~ vV/ Lp') terms in the Navier- 
Stokes equation, via the kinematic viscosity of the fluid, 
z^, and its typical velocity and length scales, V and L 
[3l 1^^ . Phenomenologically superfluids are described as a 
mixture of the superfluid and normal components with 
separate velocities, v^ and v^^ densities ps and pn, and 
viscosities z^s = and v^ > 0. The equation for Vs{r^t) 
is similar to the inviscid Euler equation [5], 



r.h.s., which describe the dissipative (ex a) and reactive 
(ex Of') mutual friction between the normal and super- 
fluid components, mediated by the superfluid vorticity 
u?s = V X Vg. Here \i is the chemical potential and 
Co = cjs/cjg. Setting Vn = 0, the ratio of the inertial 
and dissipative terms in Eq. ^ can be understood as 
the superfluid Reynolds number [6 , 



Rca (T) = [1 - a\T)]/a{T) ^ l/a{T) , 



(2) 



dt 
but with additional terms (ex v^ = v. 



a'vsXuJs^acJxlLJsXVs], (1) 
Vn) on the 



which is independent of velocity, determined only by the 
temperature dependence of a and a', which rapidly de- 
crease towards low temperatures. In the limit T ^ 0, 
Rca -^ l/o^ ^ 1, and the superfluid dynamics is ex- 
pected to become turbulent. 

This is really the case in superfluid ^He, where tur- 
bulence at the lowest temperatures is ubiquitous in mea- 
surements with moving vortices, including so-called spin- 
up or spin- down measurements, where a sudden change 
is applied to the angular velocity Q{t) of the rotating 
cryostat (see Ref. [7 for earlier experiments in a cylin- 
der or Ref. [8 in a cube). For simplicity, we restrict the 
step change to occur from to some Q or from ft to 
0. In ^He-B the viscous normal component is practically 
locked to solid-body rotation with the container, while 
the superfluid component adjusts much slower, coupled 
only by the small mutual friction from vortex motion 
with respect to the normal component. Independently of 
the shape of the container, in experimental work [71 18] 
the turbulent response of the superfluid component has 
been interpreted as the creation of a turbulent Ekman 
boundary layer adjacent to the container wall, which ex- 
pands toward the center of the container. This is similar 
to what happens in rotating viscous fluids [3]. In short, 
the transition to turbulent vortex flow in superfluid ^He 
generally resembles that in viscous fluids. 

In ^He-B the situation is different: Our non-invasive 
NMR measurements show that in a cylinder spin-down 
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FIG. 1: NMR measurement of spin down and spin up at 29 bar pressure. (Left) Calibration of the cf peak in the NMR 
absorption spectrum at 0.22 Tc. The height of the peak is plotted on the left vertical axis, its shift A/ = fd — /l from the 
Larmor frequency /l = 1.9 MHz on the right axis, and the corresponding drive Q at constant rotation in the vortex- free state 
on the bottom axis. In the insert examples of NMR spectra are shown. (Middle) Spin-down and spin-up responses plotted 
versus time. Using the calibration in the left panel, the heights and/or frequen cy shift s of the cf peak during its decay have 
been converted to corresponding Q values. The results have been fitted to Eqs. ( [4a| or (4b), to extract the characteristic time 
constants r^ and r^ in Eq. (Eq). The extrapolation of the fits to t = gives Qq for spin down or sQo for spin up. (Right) 
Spin-down and spin-up in terms of Qqt = l/{2a) versus normalized inverse temperature, Tc/T: our measurements (filled 
triangles) are compared to l/(2a) extrapolated from data on a{T) measured above 0.35 Tc in Ref. J^ (dashed line), and to 
the calculated spin down in the cylinder with 77 = 2°, which uses the mutual friction from Ref. [12] (open triangles). 



is laminar even at 0.2 Tc and Re^ ~ 10^. The same ap- 
plies for spin up in a specially prepared situation below 
0.3 Tc- In contrast, the response is turbulent when two 
quartz tuning fork sensors in the middle of the cylinder 
obstruct cylindrically symmetric flow. In the low tem- 
perature limit, where the transport of normal excitations 
is ballistic, the main difference of ^He-B from ^He is a 
much larger vortex core radius, a ~ 10...80nm. Thus 
the stability of laminar flow and the transition to tur- 
bulence depend on the quantum properties of these two 
superfluids. 

In a long cylinder, oriented parallel to the rotation 
axis, laminar rotating flow of the superfluid component 
is solid-body-like, {v^) = ft{t) x r, confined to the az- 
imuthal plane. This requires that vortex lines are highly 
polarized along the cylinder axis and their motion is 
predominantly 2-dimensional along spirally contracting 
(spin- up) or expanding (spin-down) trajectories in the 
azimuthal plane. The motion is obtained from Eq. (IT]) 
which simplifies in terms of its radial dependence to 

dn{t) 



dt 



2an{t)[no - n{t)] . 



(3) 



The solutions are for spin down ^^^(t) and spin up Q^{t): 

%(i) = no/[i + t/n], (4a) 

n^{t) = £l^o/[£ + (l-£)exp(-i/r^)], (4b) 

r = Tiy{T) = T^{T) = [2a{T) ^0]"^ (4c) 

In spin down the vortex density is oc Qq and in spin up 
oc eClo at t = 0, when the stable value of drive is reached 
after the sudden change in rotation. To conclude, laminar 
vortex flow displays a monotonic response where (vs) and 
the total vortex length L smoothly change towards the 
final state. This is to be distinguished from a turbulent 
response where L first builds up to an overshoot and then 
relaxes faster than in the laminar case. 



Measurement :-The vortex flow response is recorded 
with noninvasive NMR techniques in a smooth-walled 
quartz cylinder of 6 mm diameter. The cylinder is aligned 
along the rotation axis within < 1° [9 . The 110 mm long 
top NMR section is separated with a flat division wall 
and a small orifice from the lower part of the cylinder 
which houses two quartz tuning fork oscillators for tem- 
perature measurement [10]. Established procedures exist 
for preparing the upper NMR volume in rotating states 
with only vortex- free flow [TT]. Laminar vortex flow is 
distinguished by a shift of the NMR absorption to the so- 
called counterflow (cf) peak in the NMR spectrum [11]. 
Its height and shift from the Larmor frequency /l is a 
function of the azimuthal large-scale cf and is calibrated 
experimentally at different values of constant rotation Q 
in the vortex-free state, as seen in Fig. [T}4. 

Fig. [TjS shows an example of the decay of the cf 
peak during spin down and spin up, as expressed via 
the equivalent Q from the calibration plot in Fig. [T}4. 
During the initial sudden change of the drive Q{t) (at 
I Q 1= 0.03rad/s^) the cf peak grows continuously and 
reaches maximum height at t = when the final constant 
value of drive is attained. From there on the peak de- 
creases smoothly and monotonicly. If the vortices would 
lose their polarization because of twisting on the global 
scale, then flow along the cylinder axis would be created, 
which would lead to the formation of a new peak close to 
/l- However, except for the cf peak from azimuthal flow, 
no other peaks from large-scale flow have been observed. 
If a turbulent vortex tangle would fill the cylinder, with a 
corresponding large loss in vortex polarization and global 
flow, then all sharp cf peaks would disappear and a broad 
absorption increase close to /l would be detected. Note 
that the integrated absorption in the NMR spectrum is 
a constant at a given temperature. Thus a shift of ab- 
sorption from the cf peak to other frequencies would be 
seen as faster decay of the azimuthal flow. The slow and 
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FIG. 2: Numerical calculations of spin-down with Qo = 0.5 rad/s. 



4.3 -10 '^, and a^ = 9.1 • 10 ^ (which corresponds 



to 0.22 Tc and 29 bar). (Left) Sphere of radius 3mm; prepared to be initially in the equilibrium vortex state with A/'(0) = 
A/'eq(r2o) = 325 vortices; total number of vortex reconnections in the bulk A/'b(oo) ^ 740 and with the container surface 
Ns{oo) ~ 4. The normalized curves in the main panel represent: a) total vortex length, L(t)/L(0) - solid (blue) line; b) mean 
vortex polarization along the symmetry axis of the container, (p^) ~ dotted (green) line; c) cumulative number of the vortex 
reconnections in the bulk Nh{t)/[Nh{oo) + A^s(oo)] and on the surface Ns{t)/[Nh{oo) + A^s(oo)] - dashed (red) lines. Inserts: 
top views at t = and 100 s; top and side views at t = 1000 s. (Middle) Cylinder with diameter = height = 6 mm; inclined by 
ry = 2° from the rotation axis; A/'(0) = A/'eq(Qo) = 413, A^b(oo) ^ 10700, and A^s(oo) ^ 860. (Right) Same cylinder at different 
inclination angles 77. The time dependence of L(t)/L(0) is plotted on logarithmic scales. The overall flow behavior is laminar 
for T] < 30° and turbulent for 77 > 40°. Insert: total number of all reconnections as a function of the inclination 77. 



continuous cf absorption response can only result from 
smoothly decaying azimuthal global flow and proves the 
laminar nature of the vortex flow. 

Fig. [ip shows the characteristic decay time of lam- 
inar vortex flow as a function of inverse temperature, 
as extracted from measurements in the rotation range 
0.7. ..2. 5 rad/s. The data agree well with earlier mea- 
surements [12] of dissipative mutual friction a(T), con- 
sidering that temperature measurement in the two ex- 
periments also carries uncertainties. Spin-up times are 



longer, which is not expected from Eq. (4c), but can be 



explained by the additional time required for generating 
new vortices. The spin- up model in Eq. (4b) assumes 



that there exists an instantaneous ample supply of re- 
manent vortices covering homogeneously the cylindrical 
wall. This requires that a spin-up measurement is per- 
formed soon after a previous spin-down run, or else the 
remanent vortices are reduced to a few and spin-up starts 
from a localized vortex formation event [JJ^, followed by 
axial vortex motion with a propagating front [9]. 

To conclude, the response times in Fig. [ip constitute 
the first observations of laminar spin up and spin down, 
which remain stable with no additional dissipation be- 
yond the mutual friction from radial vortex motion up to 
large Re^ and possibly even to the limit Re^ -^ 00. 

Numerical simulation:-To shed more light on the 
laminar vortex response, we use vortex filament calcula- 
tions [11 . Since vortex creation is problematic in spin- up 
calculations, only spin down will be considered. 
• Spherical container: Fig. [2}4. The inserts show top 
and side views of calculated vortex configurations during 
laminar spin down. These configurations expand in time, 
with individual vortices moving along spiral trajectories 
to the equatorial wall, where they finally annihilate as 



half rings with a radius comparable to the core diame- 
ter. During all evolution the configuration is smoothly 
laminar, with no sign of tangling. This case resembles a 
decelerating rotating neutron star J3] and represents a 
most symmetric flow environment. 
• Circular cylinder: Fig. [2]S. To mimic a small pertur- 
bation from axial symmetry, we tilted the rotation axis 
by 77 = 2° from the cylinder axis [14 . Thus initially 
at t = the vortices are slightly inclined with respect 
to the cylinder. This aberration is small, the vortices 
are almost rectilinear, of equal length, and their density 
is constant, 2Q{0)/Hi^ as required for solid-body rotation 
(see top view in the left-most insert; hi = h/{2m^) is the 
circulation quantum). The insert in the middle shows the 
top view a little later, when vortices in the outermost ring 
have partly reconnected at the cylindrical wall, are of re- 
duced length, and are gradually annihilating. Meanwhile 
vortices inside the outer ring have changed only little in 
shape. Much later (top and side views, right-most in- 
serts) the remaining vortices are still evenly distributed 
and smooth on the scale of the inter- vortex distance L 

In the main panel of Fig. [2]S the total vortex length 
L(t) accurately follows the spin-down dependence of 
Eq. (4a), with r^ shown in Fig. [It in good agreement 
with our measurements. The mean polarization remains 
continuously high, Pz{'t) > 0.95, which applies for all tem- 
peratures > 0.20 Tc, but in particular at 0.2... O.ST^ Pz{t) 
does not change with temperature. Notice that the total 
number of reconnections on the cylindrical wall, 7Vs(oo), 
is roughly twice larger than the initial number of vor- 
tices, Af{0). This is much more than in the case of the 
sphere and suggests that the main mechanism facilitating 
vortex annihilation (in the tilted cylinder) is the surface 
reconnection on the cylindrical wall. Reconnections be- 



tween vortices in the bulk are 10 times more abundant, 
but are concentrated to the outermost ring of vortices at 
early times (where the vortex density is at maximum), 
as seen in the insert in the middle (at t = 100 s). Thus 
this numerical example is consistent with our measure- 
ments and the claim that laminar spin-down is stable in 
the bulk volume with respect to finite-size perturbation 
even at large i^e^, while reconnections and additional 
dissipation from the annihilation is concentrated in the 
outermost ring of vortices next to the cylindrical wall. 
• Tilted cylinder: Fig. |2p. By varying the inclination 
7^, we introduce a controllable perturbation. It forces the 
superfiuid to adjust to the elliptically distorted boundary 
condition and to edge effects at the tilted top and bottom 
surfaces of the cylinder. This makes it possible to identify 
a transition from laminar to turbulent spin down from the 
time dependence of L{t) at different r] in the main panel. 
When T] < 30°, L{t) displays no or a minor overshoot, 
and is followed by a slow fall-off with t~^ dependence, 
typical for laminar fiow in Eq. (4a). When r] > 40°, the 



overshoot grows larger and the spin-down decay falls on 
the fast t~^/^ dependence, typical for turbulent fiow [15 . 

The overshoot in L(t) arises when vortices reconnect 
in the bulk volume and form a turbulent tangle. Here the 
polarization along the rotation axis is largely lost and the 
kinetic energy of the global azimuthal fiow is converted 
to an increase in vortex density via frequent reconnec- 
tions. In the insert of Fig. |2p the cumulative number 
of reconnections in the bulk volume is seen to jump up 
by an order in magnitude, when rj > 40°, and simul- 
taneously the mean curvature radius of the vortices de- 
creases abruptly. The overshoot is followed by rapid tur- 
bulent decay, displayed by the overall evolution in L{t)^ 
although the turbulent tangle may not fill the cylinder 
homogeneously, but concentrates to its outer regions. It 
is surprising that such a large inclination as 7^ ~ 30... 40° 
is required before the overall spin-down response becomes 
turbulent. This observation underlines the robustness of 
laminar fiow in a cylindrically symmetric environment. 
• Cubic container: Similar calculations on a cube show 
that the spin down response is already turbulent at a low 
inclination of 77 = 2° of the cube faces with respect to the 
rotation axis [16 . The homogeneity and isotropy of the 
turbulence increase with increasing tilt angle, suggesting 
that homogeneous turbulence could become possible. 

Kelvin wave excitations:-In Fig. [2]S one can see 
bending of vortex lines which can be interpreted as 
Kelvin waves (KWs). To clarify the role of KWs dur- 
ing spin-down in the weakly tilted cylinder, compare the 
kinetic energy (per unit mass) of laminar rotation £^hd = 
{nR)^/4 with the KW energy E^^ c^ Ai^n{l - (p^))/2, 
where A = ln(^/a) c:^ 10 in ^He-B. One finds 

i?Kw/HD = E^^/E^^ ^ 4A(1 - {pI))/M. (5) 

The polarization is (p^) c^ 0.9 and the number of vor- 
tices J\f ~ 400... 40 during the evolution shown in Fig.l2]S. 



Thus the ratio i?Kw/HD varies between 0.015 and 0.15 and 
is not changing with temperature, according to the sim- 
ulation results. This means that KWs practically do not 
contribute to the total energy of the system. A second 
argument is derived comparing dissipation of KWs via 
mutual friction and via the nonlinear L'vov-Nazarenko 
KW-energy cascade [17 . One finds that the cascade may 
develop only if aK <C (1 — {pl)Y^ which is not the case 
here. We thus conclude that KWs also play no role in 
the total energy loss, which is consistent with the exper- 
imental result. 

Conclusions :-In an oriented cylinder vortex fiow re- 
sponses are laminar in the absence of strong fiow pertur- 
bations and surface pinning. This fact is deduced from 
the NMR signature of the continuously and monotonicly 
decaying superfiuid counterfiow which circulates in the 
azimuthal plane of the cylinder. The characteristic de- 
cay time has the exponential temperature dependence of 
the ^He-B mutual friction dissipation and is exclusively 
accounted for by radial vortex motion. The distinctive 
feature of such laminar fiow is a low incidence of vor- 
tex reconnections and of Kelvin wave excitations in the 
bulk volume. As a result, at 0.2 Tc the dissipation is two 
orders of magnitude lower than in the axially inhomo- 
geneous spin up in the form of a propagating turbulent 
vortex front [9 . Moreover, in the laminar case dissipa- 
tion vanishes exponentially in the limit T ^ 0, while in 
the turbulent case it approaches a finite constant value, 
as confirmed in measurements on the moving front [9] 
and on the decay of turbulent vortex tangles [18 . 
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